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The recent advances in creating nearly degenerate quantum dipolar gases in optical lattices are 
opening the doors for the exploration of equilibrium physics of quantum systems with anisotropic 
and long-range dipolar interactions. In this paper we study the zero- and finite-temperature phase 
diagrams of a system of hard-core dipolar bosons at half-filling, trapped in a two-dimensional optical 
lattice. The dipoles are aligned parallel to one another and tilted out of the optical lattice plane by 





means of an external electric field. At zero-temperature, the system is a superfluid at all tilt angles 
0 provided that the strength of dipolar interaction is below a critical value V.(@). Upon increasing 
the interaction strength while keeping @ fixed, the superfluid phase is destabilized in favor of a 
checkerboard or a stripe solid depending on the tilt angle. We explore the nature of the phase 
transition between the two solid phases and find evidence of a micro-emulsion phase, following the 
Spivak-Kivelson scenario, separating these two solid phases. Additionally, we study the stability of 
these quantum phases against thermal fluctuations and find that the stripe solid is the most robust, 
making it the best candidate for experimental observation. 


I. INTRODUCTION 


Experimental progress in trapping and controlling ul- 
tra cold atoms and molecules has led to the observation 
of magnetic and electric dipolar interactions in a vari- 
ety of systems [1-13]. These systems can be used as 
quantum simulators to study quantum transport and dy- 
namical and equilibrium properties of models featuring 
long-range and anisotropic dipolar interactions. Long- 
range dipolar interactions have been predicted to sta- 
bilize a plethora of exotic quantum phases such as p- 
wave superfluids, superfluids of multimers, charge density 
waves, stripe solids, and supersolids. Moreover, dipolar 
interactions play an important role in many models of 
strongly correlated systems, excitons with spatially sep- 
arated electrons and holes [14-16], and frustrated quan- 
tum magnets [17-19]. 

The recent success in creating a gas of polar molecules 
in an optical lattice is the first step towards the realiza- 
tion of a low-entropy initial state which is a key ingre- 
dient for any quantum simulator [7, 20]. The ability to 
change the alignment of the dipole moment by applying 
an external field has sparked numerous studies on the 
many-body phases of dipolar gases. Many-body bosonic 
and fermionic dipoles with dipoles aligned perpendicular 
to the lattice or parallel to the lattice have been studied 
extensively [21-29]. However, the theoretical and numer- 
ical effort to study dipolar gases with dipoles aligned at 
arbitrary tilt angles has either focused on continuous sys- 
tems [30-36] or has employed non-exact methods such as 
functional renormalization group [37], mean field theory 
and variational approaches [38-40]. 

In the following we present a systematic study of a 
system of hard-core, dipolar bosons trapped in a two- 
dimensional lattice. We treat the alignment of the dipoles 
as a parameter which can be adjusted via the application 
of an external field. We use the worm algorithm [41], a 





type of Path Integral Quantum Monte Carlo method, 
and study the zero- and finite-temperature phase dia- 
erams of the system at half-filling as a function of the 
tilt angle and the strength of dipolar interaction. We 
find that at low interaction strength the system is in the 
superfluid state for any value of the tilt angle. Upon in- 
creasing the interaction strength the superfluid phase is 
destabilized in favor of either a checkerboard or a stripe 
solid depending on the tilt angle. We find evidence of 
the Spivak-Kivelson scenario [23, 42—44] in the transition 
region between the two solid phases where an emulsion 
phase of the checkerboard and stripe solids is observed. 
Additionally, we study the robustness of these quantum 
phases against thermal fluctuations and show that the 
solid phases survive at temperatures higher than the crit- 
ical temperature for the disappearance of the superfluid 
phase. In particular, due to the anisotropy of the dipolar 
interaction, the stripe solid turns out to be the most ro- 
bust phase, making it the best candidate for experimental 
observation. We give predictions for actual experimental 
setups and temperatures required to observe solid phases. 
This paper is structured as follows: in section I we 
discuss the Hamiltonian describing the system. In sec- 
tion II we present the zero and finite-temperature phase 
diagrams. In section III we discuss how harmonic con- 
finement affects our results. In section IV we explore 
possible experimental realizations and provide tempera- 
ture estimates. Finally, section V concludes the paper. 


Il, HAMILTONIAN 


We study a system of hardcore, dipolar bosons with 
induced dipole moment d, confined by a two-dimensional 
optical lattice with lattice constant a and by an external 
harmonic trap. A schematic representation of this setup, 
in the absence of the external confinement, is shown in 


figure 1. In this system, the dipoles are aligned parallel 
to each other along the direction of polarization which 
forms a tilt angle 9 with an axis perpendicular to the 
plane of the optical lattice. The system is described by 
the Hamiltonian 





FIG. 1. (Color online) Schematic representation of the sys- 
tem. Dipoles are trapped in a two-dimensional optical lattice 
with their dipole moments aligned with the external electric 
field. The dipole moments make an angle 6 with respect to 
the axis perpendicular to the plane of the optical lattice. The 
anisotropic dipolar interaction depends on the angle a;; be- 
tween the direction of polarization and the relative position 
of particles. 





H=-JY~ alas t+ VS> (1 — 3008? ayy), 
r?. 
(i9) i<j 


= do Himi. (1) 


where the first term describes the kinetic energy of the 
system with hopping energy J and the second term is the 
dipole-dipole interaction with strength V = d?/a® and 
rij = |; —7;|. ayy characterizes the angle between 
the direction of the polarization and the relative position 
of the two particles given by 7;;. Here, al (a;) are the 
bosonic creation (annihilation) operators with the usual 
commutation relations and al? =, n= aay. We 
use (...) to denote nearest neighboring sites. Finally, 
= By pam weé?, where we and €; are the strength of 
harmonic confinement and the coordinate of site 2 along 
axis €, respectively, and p is the chemical potential which 
sets the total number of particles. 











Ill. ZERO- AND FINITE-TEMPERATURE 
PHASE DIAGRAMS 


In this section, we present the zero- and finite- 
temperature phase diagrams of the system described by 
Eq. 1 at half-filling and in the absence of an external 
harmonic confinement. The chemical potential is set to 
ensure the half filling condition, n = N/Ngites = 0.5. 
The effect of the harmonic confinement is discussed in 
section III. Our results are based on Path Integral Quan- 
tum Monte Carlo (QMC) simulations using the worm al- 
gorithm (WA). Unless otherwise noted, we simulate sys- 
tems with spacial dimensions L x L, where L = 24, 30, 36, 
and 42, with Neites = L x L. To extract the ground state 
phase diagram, we work at inverse temperature 6 = L/J 
which ensures that the system is effectively at zero tem- 
perature . 

The main panel in Figure 2 shows the zero-temperature 
phase diagram in the V/J — @ plane which features three 
phases: a superfluid phase (SF), a checkerboard solid 
(CB) phase, and a stripe solid (SS) phase. The superfluid 
phase possesses off-diagonal long-range order and is char- 
acterized by finite superfluid stiffness p,, which can be ex- 
tracted from simulations by measuring the winding num- 
ber in space ps = (W7?) /dL*~26 [45]. The diagonal order 
in the solid phases is characterized by a finite value of 
structure factor S(k) = >), .,explik(r —r’)| (mene) /N, 
where k is the reciprocal lattice vector. We use k = (7,7) 
and k = (0,77) to identify the CB and SS phases, respec- 
tively. 

At low interaction strength, the system is in a super- 
fluid phase for any value of the tilt angle 6. The super- 
fluid phase is destabilized towards a solid phase as the 
interaction is increased above a critical value V.(@)/J. 
Filled squares mark the onset of the CB solid, which 
forms at lower 0, while filled triangles mark the onset 
of the stripe phase which appears at larger @. 

Using energy considerations, it is easy to see that the 
anisotropic nature of the dipolar interaction leads to the 
stabilization of the CB solid at 6 = 0 and the SS phase 
at 0 = 1/2. Moreover we observe that at larger values of 
the tilt angle the superfluid phase is less stable against in- 
creasing the dipolar interaction strength. Indeed, at a tilt 
angle 6 = 7/2, the dipolar interaction strength needed to 
destroy superfluidity in favor of a solid is about a factor 
of two smaller than what needed at a tilt angle 0 = 0. 
At intermediate tilt angles, 6 ~ 7/6, the competition be- 
tween the two solid orders renders both unstable. This 
adds to the robustness of the superfluid phase, where a 
dipolar interaction V/J ~ 6 is required in order to de- 
stroy the off-diagonal long range order. 

We find the solid-SF transition points by studying fi- 
nite size effects on p, and S(k). We determine the error 
bars for the solid to superfluid transition by monitor- 
ing the order parameters characterizing the two phases, 
namely p, and S(k). We assure that the order parame- 
ter characterizing the solid phase vanishes with increas- 
ing system size as we cross the phase boundary to the 








superfluid phase and vice versa. 
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FIG. 2. (Color online) Zero-temperature phase diagram of 


the system described by Eq. 1 as a function of tilt angle 0 
and V/J, and in the absence of harmonic confinement. The 
system features a checkerboard solid (CB), a stripe solid (SS), 
and a superfluid phase (SF). The CB-SS transition follows 
a Spivak-Kievelson scenario where stripes of one phase are 
surrounded by the other phase. The shaded region indicates 
the area where we have observed this micro-emulsion phase. 


Given the nature of the interaction and the dimension- 
ality of the system, a first order phase transition between 
two ordered states is not be allowed [23, 42-44] and a 
second order phase transition is highly unlikely. First we 
investigate the nature of the phase transition between 
the CB and SS phases V/J = 8. To this end we per- 
form scans along the cut marked I on the main panel 
of figure 2. We find that the CB-SS phase transition at 
T = 0 is neither first-order nor second order and follows 
the Spivak-Kivelson micro-emulsion scenario as predicted 
in [23, 42-44] which is indicated by the shaded region in 
Fig. 2. We provide examples of the density map of the 
system in this region in figure 3(a). Here, each circle cor- 
responds to a different site, and its radius is proportional 
to the local density. In this region the system features 
a micro-emulsion phase, where stripes of CB phase are 
embedded in the SS phase or vice versa. 

We have also investigated the nature of the CB-SF 
phase transition by performing a similar scan along 
V/J = 4 line labeled as II, for a system with L = 100. 
In figure 3(b) we plot the structure factor S(7,7) as a 
function of # using filled circles. The superfluid density 
Os; is shown using filled triangles. We observe hysteresis 
in both S(z,7) and p, with a width of the hysteresis loop 
of only A@ ~ 0.015. While this observation suggests that 
the system undergoes a weak first-order phase transition, 
this is not consistent with theoretical predictions stated 
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FIG. 3. (Color online) The behavior of the system along cuts 
I and II. (a) The density distribution of the system along cut 
I through the shaded region in figure 2 at two values of 6 
shows evidence of micro emulsion phase. Here, each circle 
corresponds to a different site, and its radius is proportional 
to the local density. (b) The CB structure factor S(z,7) and 
superfluid density p,; as a function of 6 along cut II. We find 
evidence of hysteresis which indicates that this transition be- 
haves as a weak first-order phase transition for finite systems. 


above. Moreover our inability to observe a micro emul- 
sion phase intervening between the superfluid and the CB 
phase for system sizes up to L = 100 suggests that we 
may be far away from the thermodynamic limit. Hence 
if the micro-emulsion phase exists, the size of emulsion 
would be larger than this lattice size, and correspond- 
ingly large systems should be used in order to observe it 
experimentally. 


Next we present an investigation of the robustness 
of the quantum phases described above against thermal 
fluctuations. As expected, we find the solid phases to be 
the most robust against thermal fluctuations. We have 
performed scans over V/J at two tilt angles 6 = 0.271 
and 6 = 0.84 corresponding to SF-CB and SF-SS transi- 
tions at zero temperature. Our results for 6 = 0.271 are 
summarized in figure 4. ‘The phase boundaries are ex- 
tracted from finite size scaling analysis of p, and S(z,7). 
We were not able to resolve the phase boundaries within 
the shaded region with the system sizes considered in 
this paper. Within this region we expect the system to 
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FIG. 4. (Color online)(a) Finite temperature phase diagram 
at tilt angle 0 = 0.271. Upon increasing the temperature, 
thermal fluctuations destroy the CB and SF order in favor of 
a normal fluid. The CB phase melts via a two-dimensional 
Ising transition. The inset shows the scaled structure factor 
with 26/v = 0.25 for LD = 18, 24, 36, 42, 48 along cut II, at 
V/J = 4. The crossing determines the critical temperature 
T./J = 0.812 + 0.002. The SF-normal fluid transition is a 
Kosterlitz-Thouless phase transition. In (b) we show ps as a 
function of T/J for L = 18, 24, 36, 42, 48, and 66 along cut 
I at V/J = 2. The dashed line is given by T’'/z. In (c) we use 
the intersection points between the T’/7 line and the pz vs. 
T/J curves at each L to extract T./J ~% 0.66. 


undergo either a weak first order phase transition or fea- 
ture a micro-emulsion phase at zero temperature. Upon 
increasing the temperature, thermal fluctuations destroy 
the SF phase in favor of a normal fluid via a Kosterlitz- 
Thouless (KT) transition [46]. In figure 4(b) we show 
ps as a function of T/J for L = 18, 24, 36, 42, 48, and 
66 at V/J = 2, indicated as cut I in figure 4(a). In the 
thermodynamic limit, a universal jump is observed at a 
the critical temperature given by ps = 2mkgT./th?. In 
a finite size system this jump is smeared out as seen in 
figure 4(b). We extract the critical temperature using the 
finite size scaling procedure described in [47]. The dashed 
line in figure 4(b) corresponds to p, = T'/7 and its in- 
tersection points with each p, vs. T’/J curve are used to 
find T,, as shown in figure 4(c). We find T../J = 0.66 at 
V/J = 2. The CB solid melts in favor of a normal fluid 
via a two-dimensional Ising transition. We use standard 
finite size scaling as shown in the inset of figure 4 where 
we plot the scaled structure factor S(z,7)L?°/", with 
23/v = 0.25 as a function of T/J for L = 18, 24, 36, 42, 
48 along cut II, at V/J = 4. The crossing indicates a 
critical temperature T,/J = 0.812 + 0.002. 

We have performed a similar analysis at fixed tilt an- 
gle 9 = 0.84 where the zero-temperature phase diagram 
features the SF and SS phases. The results are shown 
in figure 5(a). Critical temperatures are found using the 
same methods as described for the SF and CB phases. It 
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is worth noting that while both the SS and CB phases 
are more robust against thermal fluctuations compared 
to the SF phase, at any given V/J the SS phase has crit- 
ical temperature roughly twice that of the CB phase. 

Finally, we look at the finite temperature properties 
of the CB-SS micro-emulsion phase at V/J = 8. In fig- 
ure 5(b) we show S(k) for the CB and SS phases as a 
function of T’/J at 0 = 0.48 and L = 100. We use squares 
and circles to show S(z,7) and S(z,0), respectively. The 
filled and open symbols correspond to two different ini- 
tial conditions chosen for the simulations at each T’/J. 
Since the equilibrium density distribution of the micro- 
emulsion phase is affected by the choice of initial condi- 
tions we observe large fluctuations in S(k). These results 
are in agreement with what reported in [23]. We find that 
at T’/J = 1.2 the system is in the normal fluid phase, with 
both SS and CB phases disappearing. If T/J < 1.2 the 
micro-emulsion phase can be dominated by either solid 
order, depending on the initial conditions. 
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FIG. 5. (Color online)(a) Finite temperature phase diagram 
at tilt angle 6 = 0.84. Critical temperatures are found using 
the same methods as described for the data in figure 4. (b) 
Finite temperature properties of the CB-SS micro-emulsion 
phase at V/J = 8, 0 & 0.48, and L = 100. Squares and 
circles correspond to S(z,7) and S(z,0), respectively. Filled 
and empty symbols refer to different initial conditions. At 
T/J = 1.2 the system is in the normal fluid phase, while at 
T/J < 1.2, the micro-emulsion phase is dominated by either 
solid order, depending on the initial conditions. 


IV. HARMONIC CONFINEMENT 


In a typical experimental system particles are subject 
to the optical lattice potential as well as an external har- 
monic confinement. The effect of the harmonic confine- 
ment can be taken into account in the QMC simulations 
through a site-dependent chemical potential as shown in 
Hamiltonian (1). 

Experimentally, using a sufficiently deep optical lattice 
suppresses double occupancy. In this limit, the hard-core 
condition is valid for all tilt angles. The variation of the 
chemical potential provides a scan over density hence re- 
sulting in coexistence of phases realized at half-filling (as 
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FIG. 6. (Color online) Left panels: Density distribution of 
dipoles in the presence of a harmonic potential (see text). 
Each circle corresponds to a different site, and its radius is 
proportional to the local density. Right panels: Statistics of 
the off-diagonal correlator (a! (7;)a(7;)) (see text for detail). 
A brighter color corresponds to regions of the trap where off- 
diagonal long-range order is present. ‘Top left and right pan- 
els: V/J = 2, 0 = 1.34, N ~ 511. A Mott insulator phase 
with filling factor one is observed at the center of the trap, 
while a stripe solid is stabilized in most of the outer shell. 
Bottom left and right panels: V/J = 6, 06 = 0.27, N ~ 1217. 
A superfluid phase is observed at the center of the trap, fol- 
lowed by a shell of CB solid which gives way to a superfluid 
outer shell once the harmonic potential is strong enough de- 
stroy the solid order. 





described above) with phases stabilized at other fillings. 
We use a lattice depth Vo = 50 R, where Ep is the re- 
coil energy of the molecule and a harmonic confinement 
of We, = 275 Hz. Our parameters reflect typical ex- 
perimental choices. The top and bottom left panels in 
figure 6 present the equilibrium density distribution of 
the dipoles, where each circle corresponds to a differ- 
ent site, and its radius is proportional to the local den- 
sity. In the presence of the harmonic trapping poten- 
tial we cannot use the usual winding number relations 
to measure the superfluid density in the different regions 
of the trap. However we can use the off-diagonal corre- 
lator (al (7;)a(7;)), to map the superfluid regions within 
the trap. The panels on the right-hand side of figure 6 
show the regions in the trap where the correlator is long- 
ranged. Each panel is normalized such that the brighter 
points are those exhibiting long-ranged correlations and 
thus robust superfluidity. 

The top left and right panels correspond to V/J = 2, 
tilt angle 6 = 1.34 and N ~ 511. The particles at the 
center of the trap form a Mott insulator (MI) with unit 





filling while a stripe solid is stabilized in most of the outer 
shell. The dark regions in the top right panel indicate 
an absence of off-diagonal long-range order, correspond- 
ing to the SS and MI phases. Given the steepness of 
the harmonic potential, we do not observe a superfluid 
shell separating the two insulating phases, however off- 
diagonal long-range order is present adjacent to the MI 
shell. This is due to the mismatch between the symme- 
try of the SS and the trap, which destabilizes the SS and 
allows for the build up of phase coherence in this region. 

The bottom left and right panels of figure 6 corre- 
sponds to V/J = 6 and tilt angle 6 = 0.27 for which 
the CB phase is stabilized at half filling. The total num- 
ber of particles in the trap is N ~ 1217. We observe the 
presence of superfluid phase at the center of the trap, 
followed by a shell of CB solid which gives way to a su- 
perfluid outer shell once the harmonic potential is strong 
enough destroy the solid order. The superfluid regions 
are characterized by local density different than 1 and 
0.5 as the left panel shows, while off-diagonal long-range 
order is present as seen in the right panel. On the other 
hand, the left panel clearly shows a CB structure at den- 
sity 0.5 in the inner shell with corresponding absence of 
off-diagonal long-range order (dark color in right panel). 





V. EXPERIMENTAL REALIZATION 


The system described in section | can be experimen- 
tally realized using bosonic polar molecules trapped in a 
two-dimensional lattice, provided that the lattice is deep 
enough to suppress double occupancy hence ensuring the 
hard-core condition. Once ground-state molecules are 
formed in the lattice their dipole moments can be aligned 
through the application of an external static electric field, 
realizing an interaction term of the form described in 
Eq. (1). While equilibrium states such as the ones de- 
scribed above have not been prepared yet, rapid exper- 
imental advances will allow for the realization of these 
states in the future. Our calculations suggest that, un- 
less dipoles are tilted at intermediate angles 6 ~ 7/6, for 
which superfluidity survives even for suppressed hopping 
strength J/V ~ 0.17, solid phases would be the best can- 
didate for experimental observation since they are more 
robust against thermal fluctuations. For example, for 
a system of °°K°’Rb polar molecules with permanent 
dipole moment d = 0.57D, trapped in an optical lattice 
with lattice constant a = 532 nm and depth Vo = 40E,., 
with tunneling rate J/h = 104 Hz [7], one can achieve 
J/V ~ 0.3. Hence, observation of the stripe solid would 
be possible at temperatures ~ 10 nK. With the same 
experimental setup the emulsion phase can be observed 
at temperatures ~ 6 nK. If we consider 72Na*°K with 
permanent dipole moment d = 2.7D [48] and same trap- 
ping parameters as above, one can achieve J/V ~ 0.0136 
which allow for observation of the CB solid phase and 
the SS phase at temperatures ~ 109 nK and ~ 250 nk 
respectively. It should be noted that similar behavior is 








expected in experimental systems featuring atoms and 
molecules with magnetic dipole moments, such as Ero, 
Cr, and Dy [1, 9, 13]. 


VI. CONCLUSIONS 


In summary, we presented the zero- and _finite- 
temperature phase diagrams of a system of hard-core, 
dipolar lattice bosons at half-filling as a function of the 
alignment of the dipole moments, characterized by the 
tilt-angle 6, and the strength of the dipolar interaction. 
At zero-temperature, the system features three phases: 
superfluid, checkerboard solid and stripe solid. The su- 
perfluid phase is present at all tilt angles, provided that 
the interaction strength is below V.(@)/J, and upon in- 
creasing the interaction strength the system enters one of 
the two solid phases depending on the value of 6. We ob- 


[1] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 
T. Pfau, Phys. Rev. Lett. 94, 160401 (2005). 

[2] T. Lahaye, T. Koch, B. Frohlich, M. Fattori, J. Metz, 
A. Griesmaier, 5. Giovanazzi, and T. Pfau, Nature 448, 
672 (2007). 

[3] kK. K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe’er, 
B. Neyenhuis, J. J. Zirbel, 5. Kotochigova, P. 5. Julienne, 
D.S. Jin, and J. Ye, Science 322, 231 (2008). 

[4] S. Ospelkaus, A. Pe’er, K. K. Ni, J. J. Zirbel, B. Neyen- 
huis, S. Kotochigova, P. 5S. Julienne, J. Ye, and D. S. 
Jin, Nat. Phys. 4, 622 (2008). 

[5] A. Chotia, B. Neyenhuis, S. A. Moses, B. Yan, J. P. 
Covey, M. Foss-Feig, A. M. Rey, D. S. Jin, and J. Ye, 
Phys. Rev. Lett. 108, 080405 (2012). 

[6] C.-H. Wu, J. W. Park, P. Ahmadi, S. Will, and M. W. 
Zwierlein, Phys. Rev. Lett. 109, 085301 (2012). 

[7] B. Yan, S. A. Moses, B. Gadway, J. P. Covey, K. R. A. 
Hazzard, A. M. Rey, D. S. Jin, and J. Ye, Nature 501, 
521 (2013). 

[8] T. Takekoshi, L. Reichsdllner, A. Schindewolf, J. M. Hut- 
son, C. R. Le Sueur, O. Dulieu, F. Ferlaino, R. Grimm, 
and H.-C. Nagerl, Phys. Rev. Lett. 113, 205301 (2014). 

[9] M. Lu, N. Q. Burdick, and B. L. Lev, Phys. Rev. Lett. 
108, 215301 (2012). 

[10] K. Aikawa, A. Frisch, M. Mark, S. Baier, R. Grimm, J. L. 
Bohn, D. S. Jin, G. M. Bruun, and F. Ferlaino, Phys. 
Rev. Lett. 113, 263201 (2014). 

[11] K. Aikawa, A. Frisch, M. Mark, S. Baier, R. Grimm, and 
F. Ferlaino, Phys. Rev. Lett. 112, 010404 (2014). 

[12] J. W. Park, S. A. Will, and M. W. Zwierlein, Phys. Rev. 
Lett. 114, 205302 (2015). 

[13] A. Frisch, M. Mark, K. Aikawa, S. Baier, R. Grimm, 
A. Petrov, S. Kotochigova, G. Quéméner, M. Lepers, 
O. Dulieu, and F. Ferlaino, arXiv:1504.04578 (2015). 

[14] I. V. Lerner and Y. E. Lozovik, Zh Eksp Teor Fiz 80, 763 
(1981). 

[15] A. Filinov, P. Ludwig, Y. E. Lozovik, M. Bonitz, and 
H. Stolz, J. Phys.: Conf. Ser. 35, 197 (2006). 

[16] A. K. Fedorov, I. L. Kurbakov, and Y. E. Lozovik, Phys. 
Rev. B 90, 165430 (2014). 


served signatures of a micro-emulsion phase transition, 
observing an emulsion phase of the checkerboard and 
stripe solids in the region separating the two solid phases. 
We expect solid phases to be stabilized by lattice dipolar 
fermions as well as shown in [37]. Finally, we studied 
the robustness of these phases against thermal fluctua- 
tions and showed that the solid phases can be observed 
experimentally at temperatures up to ~ 100 nk. 
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